A new light boson from MAGIC observations? 
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Recent detection of blazar 3C279 by MAGIC has confirmed previous indications by H.E.S.S. that the Universe 
is more transparent to very-high-energy gamma rays than currently thought. This circumstance can be reconciled 
with observations of nearby blazars provided that photon oscillations into a very light Axion-Like Particle occur in 
extragalactic magnetic fields. The emerging "DARMA scenario" can be tested in the near future by the satellite- 
borne Fermi LAT detector as well as by the ground-based Imaging Atmospheric Cherenkov Telescopes H.E.S.S., 
MAGIC, CANGAROO III, VERITAS and by the Extensive Air Shower arrays ARGO-YBJ and MILAGRO. 



1. MOTIVATION 

Imaging Atmospheric Cherenkov Telescopes 
(IACTs) are providing us with an impressive 
amount of information about the Universe in the 
energy interval 100 GeV — 100 TeV. Observations 
carried out by these IACTs concern gamma-ray 
sources over an extremely wide interval of dis- 
tances, ranging from the parsec scale for Galactic 
objects up to the Gigaparsec scale for the fartest 
detected blazar 3C279. This circumstance allows 
not only to infer the intrinsec properties of the 
sources, but also to probe the nature of photon 
propagation throughout cosmic distances. 

The latter fact is of paramount importance 
for very-high-energy (VHE) gamma-ray astro- 
physics, since the horizon of the observable Uni- 
verse rapidly shrinks above 100 GeV as the energy 
further increases. This is due to the fact that 
photons from distant sources scatter off back- 
ground photons permeating the Universe, thereby 
disappearing into electron-positron pairs pQ. It 
turns out that the corresponding cross section 
<t(77 — > e + e~) peaks where the VHE photon 
energy E and the background photon energy e 
are related by e ~ (500 GeV/E) eV. As far 
as observations performed by IACTs are con- 



cerned, the cosmic opacity is dominated by the in- 
teraction with ultraviolet/optical/infrared diffuse 
background photons 0, usually called Extragalac- 
tic Background Light (EBL), which is produced 
by galaxies during the whole history of the Uni- 
verse. Owing to such an absorption process, pho- 
ton propagation is controlled by the optical depth 
t(E,D), with D denoting the source distance. 
Hence, the observed photon flux <j> b s (£', D) is re- 
lated to the emitted one & em (E) by 



(1) 



Neglecting evolutionary effects on the EBL spec- 
tral energy distribution for simplicity, the optical 
depth reads t{E 1 D) ~ D/Xy(E), where Xy(E) 
is the photon mean free path for 77 — > e + e~ re- 
ferring to the present cosmic epoch. As a conse- 
quence, Eq. (JTJ simplifies as 
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The function X 7 (E) decreases like a power law 
from the Hubble radius 4.3 Gpc around 100 GeV 
to 1 Mpc around 100 TeV g]. Now, Eq. © entails 
that the observed flux is exponentially suppressed 



1 Frequency band 1.2 • 10 3 GHz - 1.2 ■ 10 6 GHz, correspond- 
ing to the wavelength range 0.25 fim — 250 fim. 
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both at high energy and at large distances, so that 
sufficiently far-away sources become hardly visi- 
ble in the VHE range and their observed spec- 
trum should anyway be much steeper than the 
emitted one. 

Yet, observations carried out by IACTs have 
failed to detect such a behaviour. A first indica- 
tion in this respect was reported by the H.E.S.S. 
collaboration in connection with the discovery 
of the two blazars H2356-309 (z = 0.165) and 
1ES1101-232 = 0.186) at E ~ ITeV [3]. 
Stronger evidence comes from the observation of 
the blazar 3C279 (z = 0.538) at E ~ 0.5 TeV by 
the MAGIC collaboration [4]. In particular, the 
signal from 3C279 collected by MAGIC in the re- 
gion E < 220 GeV has more or less the same 
statistical significance as the one in the range 220 
GeV < E < 600 GeV (6.1a in the former case, 
5.1cr in the latter). 

Turning the argument around and assuming 
standard photon propagation as described above, 
the observed spectrum ^ ^ S (E, D) can only be re- 
produced by an emission spectrum <3> cm (i?) much 
harder than for any other blazar observed so far. 

A way out of this difficulty relies upon a mod- 
ification of the emission spectrum $ em (E). A 
possibility involves the presence of strong rela- 
tivistic shocks, which can substantially harden 
$ em (E) [5] . A different option invokes photon ab- 
sorption inside the blazar, which has been shown 
to produce again an emission spectrum <E> cm (i5) 
considerably harder than previously thought [6]. 
While successful at increasing the fraction of VHE 
emitted photons, these attempts fail to explain 
why only for the most distant blazars do these 
mechanisms become important. 

A very different solution was recently proposed 
by the present authors and is usually referred to 
as the "DARMA scenario" [7 . Its characteris- 
tic feature is the presence of Axion-Like Parti- 
cles (ALPs) (more about this, later) and rests 
upon the mechanism of photon- ALP oscillation in 
cosmic magnetic fields, whose existence has def- 
initely been proved by AUGER observations [8]. 
Once ALPs are produced close enough to the 
source, they travel unimpeded throughout the 
Universe - since they do not undergo EBL ab- 
sorption - and can convert back to photons be- 



fore reaching the Earth. As a consequence, the 
effective photon mean free path A 7 . c ff (E, D) gets 
increased so that the observed photons cross a dis- 
tance in excess of A 7 (iS). Moreover, it has been 
shown that the DARMA scenario works for an 
ALP lighter than about 10" 10 eVE|. 

A deeper insight into the DARMA mechanism 
can be achieved by introducing the probability 
P 1 ^ 1 (E, D) that a photon remains a photon after 
propagation over a distance D, so that we have 



$ oha (E,D) = P^{E,D) $em(£) 



(3) 



When only photon absorption is operative, 
Eq. {2J can similarly be rewritten as 

t> ohs {E, D) = P(%(E, D) $ cm (£) , (4) 

with 

pt%(E,D)~e- D '>*W . (5) 

In the presence of photon- ALP oscillations, 
Eq. ([5]) gets replaced by 

P^(E,D) ~e- D/x ^ E) X{E,D) (6) 

and the above discussion entails X(E,D) > 1. 

Moreover, Eq. |J2J) presently becomes 

t> ohs (E, D) ~ e -n/^,ME,D) _ (?) 

with 



A T , e ff D) = 



D 



(8) 



y^-y{E,D) 

so as to guarantee consistency with Eq. (J3j) . Next 
by inserting Eq. © into Eq. §8§ we get 



X^ eS (E,D) 



D 



(9) 



Xy(E) D -X 1 (E)hxX(E,D) 

and since X(E,D) > 1 wc find X ry>e s(E,D) > 
\~(E), which is just a formal restatement of our 
previous conclusion. Still, Eq. (Q possesses the 
advantage to explicitly show that even a small in- 
crease of A 7j eff(-E, D) gives rise to a large enhance- 
ment of the observed flux <& b s (S, D). As we shall 
see, the DARMA mechanism makes \ 1>e g(E, D) 
shallower than A 7 (£'), although it remains a de- 
creasing function of E. So, the resulting ob- 
served spectrum is much harder than the one pre- 
dicted by Eq. ([2]), thereby ensuring agreement 



2 Somewhat similar ideas are discussed in [9]. 
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with observations even by adopting for far-away 
sources the same emission spectrum characteris- 
tic of nearby ones. 

Our aim is to review the main features of the 
DARMA scenario as well as its application to 
blazar 3C279. 

2. DARMA SCENARIO 

Both phenomenological and conceptual argu- 
ments lead to view the Standard Model of par- 
ticle physics as the low-energy manifestation of 
some more fundamental and richer theory of all 
elementary-particle interactions including grav- 
ity. Therefore, the lagrangian of the Stan- 
dard Model is expected to be modified by small 
terms describing interactions among known and 
new particles. Many extensions of the Standard 
Model which have attracted considerable interest 
over the last few years indeed predict the exis- 
tence of ALPs. They are spin-zero light bosons 
defined by the low-energy effective lagrangian 

I Tfl 1 ~ 

£alp = -d^ad^a-— a 2 -— F^F^ a, (10) 

where F^ v is the electromagnetic field strength, 
F^ u is its dual, a denotes the ALP field whereas 
to stands for the ALP mass d. According to 

— 1/2 

the above view, it is assumed M G F ~ 
250 GeV. On the other hand, it is supposed that 
to < G F 1/2 ~ 250 GeV. The standard Axion [TU] 
is the most well known example of ALP. As far as 
generic ALPs are concerned, the parameters M 
and to are to be regarded as independent. 

So, what really characterizes ALPs is the tri- 
lincar 7-7-a vertex described by the last term 
in £alPj whereby one ALP couples to two pho- 
tons. Owing to this vertex, ALPs can be emit- 
ted by astronomical objects of various kinds, and 
the present situation can be summarized as fol- 
lows. The negative result of the CAST exper- 
iment designed to detect ALPs emitted by the 
Sun yields the bound M > 0.86 • 10 10 GeV for 
to < 0.02 eV [TT]. Moreover, theoretical consider- 
ations concerning star cooling via ALP emission 
provide the generic bound M > 10 10 GeV, which 

3 As usual, natural Lorentz-Heaviside units with h = c = 1 
are employed throughout. 



for to < 10~ 10 eV gets replaced by the stronger 
one M > 10 11 GeV even if with a large uncer- 
tainty [10] . The same 7-7-a vertex produces an 
off-diagonal element in the mass matrix for the 
photon-ALP system in the presence of an ex- 
ternal magnetic field B. Therefore, the interac- 
tion eigenstates differ from the propagation eigen- 
states and photon-ALP oscillations show up [12] . 

We imagine that a sizeable fraction of photons 
emitted by a blazar convert into ALPs because of 
cosmic magnetic fields (CMFs), whose existence 
has been demonstrated very recently by AUGER 
observations [8]. These ALPs propagate unaf- 
fected by the EBL and we suppose that a substan- 
tial fraction of them back converts into photons 
before reaching the Earth ALPs. Owing to the 
notorious lack of information about the morphol- 
ogy of CMFs, one usually supposes that they have 
a domain- like structure [13]. That is, B ought to 
be constant over a domain of size Ldom equal to its 
coherence length, with B randomly changing its 
direction from one domain to another but keeping 
approximately the same strength. As explained 
elsewhere [H], it looks plausible to assume the 
coherence length in the range 1 Mpc — lOMpc. 
Correspondingly, the inferred strength lies in the 
range 0.3 nG- 1.0 nG [14]. 

3. PREDICTED ENERGY SPECTRUM 

Our ultimate goal consists in the evaluation 
of the probability P 7 ^ 1 (E, D) when allowance is 
made for photon-ALP oscillations as well as for 
photon absorption from the EBL. We proceed as 
follows. We first solve exactly the beam propa- 
gation equation arising from £alp over a single 
domain, assuming that the EBL is described by 
the "best- fit model" of Kneiske et al. [15] . Start- 
ing with an unpolarized photon beam, we next 
propagate it by iterating the single-domain so- 
lution as many times as the number of domains 
crossed by the beam, taking each time a random 
value for the angle between B and a fixed over- 
all fiducial direction. We repeat such a procedure 
10' 000 times and finally we average over all these 
realizations of the propagation process. 

We find that about 13% of the photons arrive 
to the Earth for E = 500 GeV, representing an 
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enhancement by a factor of about 20 with respect 
to the expected flux without DARMA mechanism 
(the comparison is made with the above "best- 
fit model"). The same calculation gives a frac- 
tion of 76% for E — 100 GeV (to be compared 
to 67% without DARMA mechanism) and a frac- 
tion of 3.4% for E = 1 TeV (to be compared to 
0.0045% without DARMA mechanism). The re- 
sulting spectrum is exhibited in Fig. 1. The solid 
line represents the prediction of the DARMA sce- 
nario for B ~ 1 nG and Ldom — 1 Mpc and the 
gray band is the envelope of the results obtained 
by independently varying B and L& om within 
a factor of 10 about such values. These con- 
clusions hold for m <C 10 _10 eV and we have 
taken for definiteness M ~ 4 ■ 10 11 GeV but we 
have cheked that practically nothing changes for 
10 11 GeV < M < 10 13 GeV. 




Photon energy (GeV) 



Figure 1. The two lowest lines give the frac- 
tion of photons surviving from 3C279 without the 
DARMA mechanism within the "best-fit model" 
of EBL (dashed line) and for the minimum 
EBL density compatible with cosmology (dashed- 
dotted line) . The solid line represents the predic- 
tion of the DARMA mechanism as explained in 
the text. 



Our prediction can be tested in the near future 
by the satellite-borne Fermi LAT detector as well 
as by the ground-based IACTs H.E.S.S., MAGIC, 



CANGAROO III, VERITAS and by the Exten- 
sive Air Shower arrays ARGO-YBJ and MILA- 
GRO. 
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